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Abstract 



We use the field theoretical model to perform relativistic calculations of 
neutrino energy losses caused by the direct Urea processes on nucleons in 
the degenerate baryon matter. By our analysis, the direct neutron decay 
in the superdense nuclear matter under beta equilibrium is open only due 
to the isovector meson fields, which create a large energy gap between pro- 
tons and neutrons in the medium. Our expression for the neutrino energy 
losses, obtained in the mean field approximation, incorporates the effects of 
nucleon recoil, parity violation, weak magnetism, and pseudoscalar interac- 
tion. For numerical testing of our formula, we use a self-consistent relativistic 
model of the multicomponent baryon matter. The relativistic emissivity of 
the direct Urea reactions is found substantially larger than predicted in the 
non-relativistic approach. We found that, due to weak magnetism effects, 
relativistic emissivities increase by approximately 40-50%, while the pseu- 
doscalar interaction only slightly suppresses the energy losses, approximately 
by 5%. 

PACS number(s): 97.60.Jd , 21.65-hf , 95.30.Cq 
Keywords: Neutron star. Neutrino radiation 
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I. INTRODUCTION 



Modern calculations based on relativistic equations of state indicate that the neutron 
star cores consist of neutrons with the admixture of protons, electrons, muons and some 
exotic particles (including hyperons, K- mesons, quarks and so on...). The composition 
is governed by the charge neutrality and equilibrium of the medium under weak processes 
Bi — >■ i?2+/+^'h B2+I Bi+i>i, where Bi and B2 are baryons (or quarks), and / is a lepton, 
either an electron or a muon. These reactions, widely known as the direct Urea processes, are 
a central point of any modern scenarios of evolution of neutron stars. Neutrino energy losses 
caused by the direct Urea processes lead to a rapid cooling of degenerate neutron star cores 
0. The corresponding reactions on nucleons, n —>■ p + I + ui, p + I n + ui, are the most 
powerful sources of neutrinos and antineutrinos in cooling neutron stars. In spite of widely 
adopted importance of these reactions, the corresponding neutrino energy losses are not 
well investigated yet. A simple formula suggested by Lattimer et al. |^ more than ten years 
ago has been derived in a non-relativistic manner. Actually, however, the superthreshold 
proton fraction, necessary for the direct Urea processes to operate in the degenerate nuclear 
matter, appears at large densities, where the Fermi momenta of participating nucleons are 
comparable with their effective mass. Moreover, according to modern numerical simulations, 
the central density of the star can be up to eight times larger than the nuclear saturation 
density [||]. This implies a substantially relativistic motion of nucleons in the superdense 
neutron star core. The appropriate equation of state for such a matter is actually derived in 
the relativistic approach, and the relevant neutrino energy losses must be consistent with the 
used relativistic equation of state. Some aspects of this problem was studied by Leinson and 
Perez |Q , who have estimated relativistic effects of baryon recoil and parity violation in the 
direct Urea processes including also effects of the baryon mass difference. This approach is 
useful, for example, for the direct Urea reactions changing the baryon strangeness because, 
in this case, the contribution of the weak magnetism into the matrix element of the beta- 
decay is not well known even for free particles. For nucleons, the relativistic regime should 
incorporate the effects of weak magnetism and pseudoscalar interaction, which drastically 
influence the neutrino emissivity of the corresponding reactions. To demonstrate this, in 
the present paper we consider the totally relativistic direct Urea process on nucleons. We 
utilize the Walecka-type relativistic model of baryon matter ^ , where the baryons interact 
via exchange of a, u, and p mesons, and perform the calculation of the neutrino energy 
losses in the mean field approximation. This approximation is widely used in the theory of 
relativistic nuclear matter, and allows to calculate in a self-consistent way the composition of 
the matter together with energies, and effective masses of the baryons. In Section II we begin 
with considering the relativistic kinematics of the neutron beta decay in the medium under 
beta equilibrium. We consider a free gas model and some field theoretical models of nuclear 
matter to demonstrate that the direct neutron decay in neutron stars is open only due to 
strong interactions caused by isovector mesons. We shortly discuss the energies and the wave 
functions of nucleons in the mean field approximation and demonstrate the nonconservation 
of the charged vector current of nucleons in this model. The matrix element of the neutron 
beta decay is derived in Section III. In Section IV we calculate the neutrino energy losses 
caused by the direct Urea on nucleons in the degenerate nuclear matter under chemical and 
thermal equilibrium. In Section V we inspect the non-relativistic limit of the neutrino energy 
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losses in order to compare this with the expressions earher obtained in and . Efficiency 
of the relativistic approach is numerically studied in Section VI. We evaluate the neutrino 
energy losses due to the direct Urea processes on nucleons in the multicomponent baryon 
matter under beta equilibrium and compare the relativistic result with that predicted by the 
known non-relativistic formula. We specially discuss the contributions of weak magnetism 
and pseudoscalar interaction. Summary and conclusion are in Section Vll. In Appendix, we 
discuss some details of the model used for the nuclear matter. 

In what follows we use the system of units h = c = 1 and the Boltzmann constant 
ks = 1- Summation over repeated Greek indexes is assumed 



II. KINEMATICS OF THE REACTION AND MODELS OF NUCLAR MATTER. 

Before calculating the neutrino energy losses caused by the direct Urea processes, let us 
examine relativistic kinematics of the reaction, n — > p+Z + z/^, in the degenerate matter under 
beta equilibrium. In what follows we consider massless neutrinos of energy and momentum 
ki = (^uJi, kij with uji = ki . The energy-momentum of the final lepton / = e~ ^jjL~ of mass 



mi is denoted as k2 = (1^2,^2) with U2 = yk2 +mf. Thus, the energy and momentum 
conservation in the beta decay is given by the following equations 

(p) - Ep {^) -uji- 002 = 

p-p-ki-k2 = 0, (1) 

where E^ (p) and Ep (j?) are the in-medium energies of the neutron and the proton respec- 
tively. 

The energy exchange in the matter goes naturally on the temperature scale ~ T, which 
is small compared to typical kinetic energies of degenerate particles. Therefore the momenta 
of in-medium fermions can be fixed at their values at Fermi surfaces, which we denote as 
Pn, Pp for the nucleons and pi for leptons respectively. Since the antineutrino energy is 



uji ~ T, and the antineutrino momentum ki 



T is much smaller than the momenta of 



other particles, we can neglect the neutrino contributions in Eqs. (|l]). Then the momentum 
conservation, Pp + pi = Pn, implies the well known "triangle" condition, 

Pp+Pl>Pn, (2) 

necessary for the Urea processes to operate. However, the momentum conservation is neces- 
sary but insufficient condition for the direct beta decay to occur. The energy conservation 
requires a one more condition 



En (pn) - Ep (pp) = ^{pn- Pp) + mf, (3) 
Squaring of both sides of this equation gives 

(En (pn) - Ep ipp)f - (pn - P^f = ml (4) 
If we denote the energy-momentum transfer from the nucleon as 
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K ={E^- Ep, pn - Pp) , 



(5) 



the Eq. @ can be readily recognized as K"^ = mf. This condition naturally arises from the 
time-like momentum of the final lepton pair, K = ki + k2 — k2, and can be satisfied only in 
the presence of the energy gap between spectrums of protons and neutrons. 

Consider now the condition of beta equilibrium, /in — /ip = /i;, where //n, /ip, and fii 
are the chemical potentials of neutrons, protons, and leptons respectively. The chemical 
potentials of degenerate particles can be approximated by their individual Fermi energies, 
yielding 

E„ (pn) - Ep (pp) = ^pf + mf (6) 

Generally^ this equation is not the same as Eq. @. However, if the "triangle" condition 
is fulfilled then the direct neutron decay is open, and one has pi = Pn — Pp- In this case the 
Eq. (I) becomes identical to the energy consrvation (|^). Thus under beta equilibrium, the 
"triangle" condition is the necessary and sufficient for the direct neutron decay to occur. It 
should be emphasized however that the stated conditions are consistent only in the presence 
of the energy gap between protons and neutrons. Therefore the possibility of the direct 
neutron decay depends substantially on the model of nuclear matter. 

A. Free gas model 

Consider, for example, a degenerate free gas consisting of neutrons, protons, and electrons 
under beta equilibrium. In this case the energy gap exists only due to the mass difference 
of a neutron and a proton. If we denote the masses as and Mp respectively, then the 

corresponding Fermi energies are (pn) = \/ Ml + p^, and Ep (pp) = + Pp. Due to 

charge neutrality the number density of electrons, oc pf, equals to the number density 
of protons, Up oc Pp. This implies Pe = Pp and the equation of chemical equilibrium 
becomes of the form 

V^e +PI + +PI = V^n + Pi (7) 

Solution of this equation 

Pp (Pn) = :\ J {pI + Ml - (Mp + m,f) {pl + Ml - (Mp - m,f) (8) 
2^Ml + pi V 

gives the proton Fermi momentum as a function of the neutron Fermi momentum in the 
beta-equilibrated gas of protons and neutrons. As required by the "triangle" condition, 
2pp > Pn, the direct neutron decay in such a medium is open when pp is larger than Pn/2. 
So, solution of the equation 



We remind that, due to the modified Urea processes, the nuclear matter attains chemical equiUb- 
rium even if the direct beta decay is forbidden. In this case the Eq. (|6|) is vaUd but pi ^ |Pn — Pp|. 
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p. iPn) = f (9) 

gives the critical value of the neutron Fermi momentum above which the direct neutron 
decay is forbidden. We find 



= ^ = = 2. 381 1 MeV. (10) 

/2M2 + 2mf - M2 



Thus the direct neutron decay is forbidden if the number density of neutrons is larger than 



n1 = = 5.932 X 10^^ cm'^ (11) 



37r2 

This number density is much smaller than that typical for neutron star cores, therefore the 
direct neutron decay in the cooling neutron stars can occur only due to strong interactions. 

Notice, some of models of strong interaction also forbid the direct neutron decay. Con- 
sider, for example, a simple model for the baryon matter 0, which contains fields for baryons 
and neutral scalar (ex) and vector (u;^) mesons. This model reproduces in a simple way the 
empirical NN scattering amplitude and the bulk nuclear properties. However, the neutral 
scalar and vector mesons equally interact with protons and neutrons, which therefore have 
identical energy spectrums. In this case, the energy and momentum in the direct neutron 
decay can not be conserved simultaneously. 

To allow the direct neutron decay, the model of nuclear matter must be generalized to 
include some additional degrees of freedom and couplings, which are able to create a large 
energy gap between possible energies of the proton and neutron. For this purpose, besides 
isoscalar mesons a and uj, the model should include also isovector mesons. The isovector 
meson couples differently to protons and neutrons, thus creating the energy gap necessary 
for the direct neutron decay. 



B. Field theoretical model. Mean field approximation. 

In the following we consider a self-consistent relativistic model of nuclear matter in which 
baryons interact via exchange of isoscalar mesons cr and u and an isovector meson p (See 
Appendix). In the mean field approximation, when the contribution of mesons reduce to 
classical condensate fields (a) = (Tq, (o;^) = Wo5'^°, (b'^'j = {0,0, po) S'^^, only the baryon 
fields must be quantized. This procedure yields the following linear Dirac equation for the 
nucleon 

{td^l^ - g^^'ujo - ^^7p7Vor3 - (M - g.ao)^ ^ (x) = 0, (12) 

Here and below we denote as Ts, and t± = (ri ± iT2) /2 the components of isospin operator, 
which act on the isobaric doublet (x) of nucleon field; M = 939 MeV is the bare nucleon 
massi. 



^As explained in previous Section, the mass difference of the proton and the neutron can be 
neglected considering the nucleon as an isobaric doublet. 
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The stationary and uniform condensate fields equally shift the effective masses 

M* = M- g^a, (13) 
but lead to different potential energies of the proton and neutron 

Un = gi^ujQ - ^gpPo, Up = g^uo + ^gpPo, (14) 

thus creating the energy gap Un — Up = —gpPo between possible energies of protons and 
neutrons. 



The exact solutions of Eq. (12) can be found separately for protons and for neutrons. 
In our case of a stationary, uniform system, solution for the neutron is a spinor plane wave 

V'n (x) = N^Un exp {-iEj + ipr) , (15) 

where the neutron energy is given by i?n (p) = \Jv^ + M*^ + t/n, and the spinor (P) should 
be found from the following equation 

((E, - f/J 7° - 7P - M*) = 0. (16) 

It is easily to understand that solution of this equation = Un (P) is a free-like spinor 
constructed from the neutron kinetic momentum 



= (E„ -U^,p) = [^Jf + M*2,pj , (17) 

For the proton, one has 

% (x) = NpUp {P') exp i-iEpt + %^r) , (18) 
where the spinor Up (P') obeys the equation 

[^p - Up) 7° - 7P - M*) Up = (19) 
and depends respectively on the proton kinetic momentum 

P"- = (Ep -UpJ) = [^^p'^ + M*^ ff^ . (20) 

In what follows we denote hj e = y/p^ + M*'^, e' = \/p''^ + M*^ the kinetic energy of the 
neutron and the proton respectively. So the normalization factors are of the form 



and the single-particle energies are E^{p) = e + U^, and Ep {ff) = e' + Up. 



At the Fermi surfaces the single-particle energies are E^ (pn) = J M*"^ + Pn + ^n, and 



Ep {pp) = ^M*2 + p'^ + Up. In this case the Eq. (|J) takes the form 

(^^Mf+pl - ^Mf+pl + U^-Upy- (p, - ppf = ml (22) 

Typically Un — Up = —gpPo ~ 100 MeV, so this equation can be readily satisfied simultane- 
ously with the momentum conservation Pn — Pp = Pi- 
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C. Nonconservation of the charged vector current of nucleons. 

The Lagrangian density (^l|) ensures a conserved isovector current [^] 

f A* = ^ij^^Tij + 6. X a^f ^ = 0, (23) 

where B^'^ = d^lf — d'^h^ is the p meson field strength tensor. Besides the directly nucleon 
contribution the conserved current (^Sf ) includes also the contribution of the isovector field 
6'*, which obeys the field equations 

+ mjfe'^ = l^^^^A'r^, d% = 0. (24) 

By the use of Eq. (p^ the condition df^T'^ = may be transformed as 

idf, (?/'7^r^) = Qpbf, x tp^T^i^- (25) 
In the mean field approximation this gives 

id^ (ij^r+^p) = -^?pPoV'7°r+^• (26) 
By introducing the covariant derivative 

(^-Wo, vj, (27) 

we can recast the Eq. (^6D to the following form 

D, {tPrT+^) = 0. (28) 
At the level of matrix elements this can be written as 

in q^^Un (P) = 0, (29) 
where is the kinetic momentum transfer 

q^ = (En - Ep + gpPo, p-p) = {s-e',p-f/) (30) 

Thus the matrix element of the transition current is orthogonal to the kinetic momentum 
transfer but not to the total momentum transfer from the nucleoni. Note that this effect 
originates not from a special form (^) of the conserved isovector current in the medium 
but is caused by the energy gap between the proton and neutron spectrums. This follows 
directly from the Dirac equation ([T^), which ensures the Eq. ( p9|) with 



q^ = P^- P'^ = (En - Ep - f/n + t/p, p - j/) , (31) 

which coincides with Eq. (^O]) because U^ — Up = —gpPo- 



^Some consequences of this fact for the weak response functions of the medium are discussed in 
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III. MATRIX ELEMENT OF THE NEUTRON BETA DECAY. 



In the lowest order in the Fermi weak couphng constant Gp, the matrix element of the 
neutron beta decay is of the form 

(/I {S - 1) |z) = -^^N^N^Ul ih) 1, (1 + 75) ^ {-h) p (P'l (0) |P), X 

X (27r)' 5 (E„ -E^-uoi- 002) 6{p-p -h- h) , (32) 

were C = cos 6c = 0. 973 is the Cabibbo factor. The effective charged weak current in the 
medium consists of the polar vector and the axial vector, J'^ (x) = (x) + {x). 

Our goal now is to derive the nucleon matrix element of the charged weak current in the 
medium. Consider first the polar-vector contribution. By the use of the isovector current 
(p3|) one can construct the conserved electromagnetic current in the medium 

JL = l^^Yi^ + + ^ (^Aa^^vi/) , a, J,^ = 0. (33) 
The last term in Eq. (|33D is the Pauli contribution, where 2a'^'^ = 'j'^'j'^ — 'j'^'j'^ and 

A = Ap^(l+r3) + A,i(l-r3). (34) 

In the mean field approximation, we replace the magnetic formfactors of the nucleon with 
anomalous magnetic moments of the proton and the neutron, Ap = 1.7928 and An = —1.9132. 

By the conserved- vector-current theory (CVC), the nucleon matrix element of the charged 
vector weak current is given by 

P (^'1 |P)„ = p (P'l J,^ |P)p - „ (P'l J,^ |P)„ . (35) 
For the electromagnetic transitions one has 

p (P'l J!f^ (0) |P)p = ^p (P') (^Y + ^Ap a^'^g.) Up (P) , (36) 



n (P'l JL (0) |P)n = (P') (^A, a^'^g.) (P) (37) 
with q^ = P — P', as given by Eq. (PDj). Thus, in the mean field approximation, we obtain 

p(P'|y^(0)|P), = np(P') 

The second term in Eq. (^), describes the weak magnetism effects. By the use of the Dirac 
Eqs. ( |T6|) and ([19|) for the nucleon spinors we find 

p(P'|g,l^'^(0)|P)^ = (39) 

in accord with Eq. (^). 
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2M 



(38) 
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Consider now the axial- vector charged current. This current is responsible for both the 
np transitions and the pion decay. In the limit of chiral symmetry, 0, the axial-vector 

current must be conserved. In the medium with p meson condensate, this implies 

lim D.AP (x) = 0, (40) 

m,r— >0 



where the covariant derivative is defined by Eq. (|27|). At the finite mass of a pion, 
m^, the axial- vector charged current is connected to the field 7r_ {x) = (tti + ^2) /a/2 of vr" 
meson. For a free space, this relation is known as the hypothesis of partial conservation of 
the axial current (PCAC). In the medium the PCAC takes the form 

D,AP (x) = mlUn^ (x) , (41) 

where m.,^ = 139 MeV is the mass of vr-meson, and /^r is the pion decay constant. 

With allowing for interactions of the pions with nucleons and p mesons the Lagrangian 
density for the pion field is of the form ^ 

^^ = \ [("^M^ - 9pK ^ ^) ■ {9^^ - 9pb^ X vf ) - mln ■ vr] + ig-^ilj-i^T ■ nip. (42) 

In the mean field approximation this results in the following equation for the field of vr" 
meson 

((2^° + g.pof - (zV)' - ml^ 7r_ (x) = -V2ig^iljpj5i^^, (43) 

where g.,^ is the pion-nucleon coupling constant. 

For the nucleon transition of our interest, the Eq. (|4l|) gives 

p(P'| q,AP (0) |P)„ = ^ml^,{P'\ n (0) |P), (44) 



Here the right-hand side can be calculated by the use of Eq. (|^). We obtain 

p(P'| q,AP (0) |P)„ = - ^f_!f u, (PO 75«n (P) . (45) 

This equation allows to derive the nucleon matrix element of the axial- vector charged current. 
Really, to construct the axial- vector matrix element of the charged current, caused by the 
nucleon transition, we have only two independent pseudovectors, consistent with invariance 
of strong interactions under T2 isospin transformation, namely: Up{P') 'j^'j^u^lP), and 
Up (P') qP'j^u^ (P). This means that the matrix element of the axial- vector charged current 
is of following general form 

p(P'| AP (0) |P)^ = Ca Up (P') {Y^, + F, g'^75) «n {P) ■ (46) 

Here, in the mean field approximation, we set Ca — 1.26, while Fq is the form-factor to be 
chosen to satisfy the Eq. (^), which now reads 

Ca (-2M* + Fq q') Up (P') 75«n (P) = - ^f^'^^t up (P') 75«n (P) • (47) 
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To obtain this we used the Dirac Eqs. ([TBI) and (^) for the nucleon spinors. Thus 

CA{2M*-F,g^) = y^^^. (48) 

'"'IT y 

In the mean field approximation, we assume that the couphng constants are independent of 
the momentum transfer. By setting = in Eq. ( ^5] ) we obtain the Goldberger - Treiman 
relation 

f^g^ = V2M*Ca. (49) 

By inserting this in (^) we find 

2M* 

F,=- T^T^y (50) 
[ml — q'^) 



Thus, with taking into account Eqs. (|38|), (^61), and (|50D, the total matrix element of the 
neutron beta decay is found to be 

G C 

Mfi = -i^=-ui {k2) 7/. (1 + 75) I' i-h) X (51) 



xu, {P') \^CyY + ^Cmct^'^I'^ + Ca (7^5 + F, g'^75) 



where, in the mean field approximation, we assume 

Cv = l, Ca/ = Ap - A„ ~ 3.7, Ca = 1.26. (52) 

Note that the matrix element obtained is of the same form as that for the neutron decay 
in a free space, but with the total momentum transfer replaced with the kinetic momentum 
transfer. Due to the difference in the neutron and proton potential energy, the kinetic 
momentum transfer 

q = P-P' = {e-e',p-p) (53) 

to be used in the matrix element (|5T1) differs from the total momentum of the final lepton 
pair 

K={e-e' + Un-U^,p-p) (54) 

This ensures > 0, while = (e — e')^ — {p — P' f' < . 

The square of the matrix element of the reaction summed over spins of initial and final 
particles is found to be§: 



"^Here we rectify an error made in the journal version of the paper. The correct expression can 
be obtained from that published in |jl^ by simple replacement Cm Cm/'^- The author is 
grateful to M. Prakash and S. Ratkovic who have pointed out this error. 
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Mfif = 32GlC^ [(Cl - C^) M*2 {k,k2) + {Ca - Cvf {k,P2) (k^Pi) 

+ iCA + Cvf{kiPi)ik2P2) 

+Cm^ [2Ca ((fciPi) {k2P2) - {kiP2) (fcaPi)) 

+Cv {{kik2) (P1P2 - - {k^P^ - kiP2) {k2Pl - k2P2))' 

'M*^kiP2){3 (fcsPs) - (fcsPi)) 

+M*2 {kiPi) (3 (A;2Pi) - (^2P2)) + (A;iA;2) (P1P2 - M*^)' 

- (A;iPi + kiP2) {k2Pi + k2P2) (P1P2)] 
+ClP, (2M* + P, (m*2 - (P1P2))) [{k,k2) (m*2 - 

- (A;iPi - A;iP2) (A;2Pi - A;2P2)]] 



'4M2 



with Pi = and P2 = {£',p')- 



(55) 



IV. NEUTRINO ENERGY LOSSES 



We consider the total energy which is emitted into neutrino and antineutrino per unit 
volume and time. Within beta equilibrium, the inverse reaction p + l ^ n + vi corresponding 
to a capture of the lepton /, gives the same emissivity as the beta decay, but in neutrinos. 
Thus, the total energy loss Q for the Urea processes is twice more than that caused by the 
beta decay. Taking this into account by Fermi's "golden" rule we have 



Q = 2/ 



d^k2d^kid^pd^p' 



7iM^rc^i/n(i-/p)(i-/o 



X 



(27r)i22cu22cui2e2e' 
(27r)' S (Pn ip) - Pp (p') -UJi-UJ2)S{p-f/-h-k2 



(56) 



Antineutrinos are assumed to be freely escaping. The distribution function of initial neutrons 
as well as blocking of final states of the proton and the lepton / are taken into account by the 
Pauli blocking-factor (1 — /p) (1 — /;). The Fermi-Dirac distribution function of leptons 
is given by 



fi M 



1 



exp {u2 - Hi) /T + 1' 
while the individual Fermi distributions of nucleons are of the form 



exp ((£ + [/„-//„) /P) + l' 



(57) 



(58) 



exp((£' + t/p-//p) /r) + i' 



(59) 



By neglecting the chemical potential of escaping neutrinos, we can write the condition of 
chemical equilibrium as iJ>i — i^n — A*p- Then by the use of the energy conservation equation. 
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£ + Un = e' + Up + UJ2 + uji, and taking the total energy of the final lepton and antineutrino 
as UJ2 + uji = fJ^i + uj' we can recast the blocking-factor as 

= /n {e) (1 -fn{e- 00')) (1 - fi (/xz + 00'- UJ,)) , (60) 

where u' ~ T. 

Furthermore, since the antineutrino energy is Ui ~ T, and the antineutrino momentum 
ki ~ T is much smaller than the momenta of other particles, we can neglect the neutrino 
contributions in the energy-momentum conserving delta-functions 

6 {e + U^- e' - Up - uji - 002) 6 (p- p -k,- ko^ 
~ 5 (e + f/n - e' - f/p - ^2) 5 (p - P - (61) 

and perform integral over cPp' to obtain p = p — k2 in the next integrals. 

Nucleons and leptons in the neutron star core are strongly degenerate, therefore the main 
contribution to the integral comes from narrow regions near the corresponding Fermi 
momenta. Thus we can set \p\ = pn, ^2 = Pi in all smooth functions under the integral. 

The energy of the final lepton is close to its Fermi energy fii = fi^ — fip. Here the chemical 
potentials of nucleons can be approximated by their individual Fermi energies = + f^n, 
fip = Ep + Up. This allows us to transform the energy-conserving (5-function as 



S{en - V^n + Pf - ^PnPl COS Oi + M*^ + U^ - Up - fn) 

^-6 (cosOi - (pl-pl+p^)] , (62) 



PnPl V 2pnP« 

where 61 is the angle between the momentum p^ of the initial neutron and the momentum 
Pi of the final lepton. Notice, when the baryon and lepton momenta are at their individual 
Fermi surfaces, the 6- function (^) does not vanish only if Pp + pi > Pn- 

Further we use the particular frame with Z-axis directed along the neutron momentum 
Pn- Then 

Pi = (0, 0, Pn, £n) 

ki = uji (sin6'i^, 0, cos 6',^, 1) 

/C2 = {pi sin 61 cos (pi, pi sin 61 sin ipi, pi cos 6^, fJ^i) (63) 
The energy-momentum of the final proton is defined by conservation laws: 

P2 = {-pi sin 01 cos ipi, -pi sin 61 sin ipi, Pn- Pi cos 6^, Sp) (64) 



Insertion of (^) and (^) in the square of the matrix element (|55D yields a rather cumbersome 
expression, which, however, is readily integrable over solid angles of the particles. 

Since we focus on the actually important case of degenerate nucleons and leptons, we 
may consider the neutrino energy losses to the lowest accuracy in T/fii. Then the remaining 
integration reduces to the factor 
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J dujiiufdiu'de /n (e) (1 -/„(£- cu')) (1 - fi (/xz + lu' - lui)) 



dijj\ 



457 



-oo exp u '/T-IJq 1 + exp (cji - w') /T 5040 
Finally the neutrino emissivity is found to be of the formi: 
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2 - £p) - (3 (£n - £p) - Pz j 



+CvCm 



2M 



M* 



+Ca ( Cv + J^Cm ) {Sn + Sp) [pf - (£n " 



+ "^16M2 



+ (p2 - (e, - spf) {2el + 25j - p^) 

- (pf - (Sn - Spf) (Sn + Spf (2^^ " 2£p - jJ^i) 



-ClM*'^ (l + m^$) - ep)' +pf)-2 - £p) pf] } 

with G (a;) = 1 if a; > and zero otherwise. In the above, the last term, with 



2 ' 



(65) 



(66) 



(67) 



represents the contribution of the pseudoscalar interaction. The "triangle" condition pp + 
pi > Pn, required by the step-function, is necessary for conservation of the total momentum 
in the reaction and exhibits the threshold dependence on the proton concentration. 



V. NON-RELATIVISTIC LIMIT 

Under beta-equilibrium, the superthreshold proton fraction in the core of neutron stars 
appears at large densities, when Fermi momenta of nucleons are of the order of their effective 
mass. Therefore the "triangle" condition is inconsistent with the non-relativistic limit. 
Nevertheless, it is necessary to show how the expression given by Lattimer et al. can be 
formally obtained from Eq. ([66|). For this purpose we consider the non-relativistic limit of 
Eq. (|66D by neglecting the "triangle" condition. 

The non-relativistic approximation is valid when pn,pp,pi -C M*. However, the smallness 
of the particle momenta is not enough to point out unambiguously the leading terms in Eq. 
. The relative contributions of various terms depend also on electron abundance in the 



^See footnote to Eq. (|5|) 
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medium. If the electron fraction is as large that M*fii ^ then we obtain the result of 
Lattimer et al. HSl: 



Qi 



4577r 
10080 



(68) 



widely used in simulations of neutron stars. However, when ~ M*fii, what actually takes 
place in the non-relativistic nucleon matter under beta-equilibrium, we obtain 



Qn 



4577r 
10080' 



(69) 



Here the additional term is due to the nucleon recoil. When the electron (proton) fraction 
is small, this contribution is comparable to the terms of Eq. (|68D. 



VI. EFFICIENCY OF THE RELATIVISTIC APPROACH 

Neutrino energy losses caused by the direct Urea on nucleons depend essentially on the 
composition of beta-stable nuclear matter. Therefore, in order to estimate the relativistic 
effects, we consider the model of nuclear matter, which includes nucleon and hyperon degrees 
of freedom (See Appendix). The parameters of the model are chosen as suggested in Ref. 
1^ to reproduce the nuclear matter equilibrium density, the binding energy per nucleon, the 
symmetry energy, the compression modulus, and the nucleon effective mass at saturation 
density uq = 0.16 fm~^. The composition of neutrino- free matter in beta equilibrium 
among nucleons, hyperons, electrons and muons is shown in Fig. 1 versus the baryon number 
density nb, in units of tt-q (left panel). The right panel represents the results of self-consistent 
calculation of the effective baryon masses and individual Fermi momenta of the baryonsi. 

At number densities smaller than the saturation density the matter consists mostly of 
neutrons with a small admixture of equal amounts of protons and electrons, however, the 
proton fraction grows along with the matter density. Near saturation density the electron 
chemical potential reaches the muon mass, and the muon fraction appears also growing along 
with the matter density together with the electron fraction up to rih ~ 0.3 fm^^. Above this 
density of nuclear matter, constituents other than neutrons, protons, electrons, and muons 
appear. The S~ hyperons appear when their lowest energy state first lies below /in + /^i- Due 
to the charge neutrality, increasing of the S~ fraction suppresses the lepton abundance. 



^In order to apply the non-relativistic Eq. ( pSj ) to relativistic nucleons the a uthors su ggest 
replacing the effective Dirac mass with the relativistic Fermi energy, M*^ — > ^Jpn,p + M*^. The 
latter is termed as "the effective Landau mass" . This trick can not be justified from the theoretical 
point of view. Firstly because the mass term and the energy term in the Dirac equation are of 
different matrix structure. Secondly, by the correct definition, the effective Landau mass has no 
relation to relativism but is intended to take into account interactions between particles near the 
Fermi surface. In this meaning the effective mass (^) is actually the effective Landau mass. 



14 



2 4 6 8 



2 4 6 8 



Fig. 1. The left panel shows individual concentrations for matter in beta equilibrium 
among nucleons, hyperons, electrons and muons as a function of the density ratio nt/nQ. The 
right panel represents the self-consistent effective baryon masses (solid lines) and individual 
Fermi momenta of the baryons (dashed lines). 

Above the density rih ~ 0.48 fm~^, beta equilibrium requires appearance of A hyperons 
because the lowest energy state for a A lies lower than /^n. At densities larger than rib — 0.91 
fm~^ fraction of S° hyperons also exists. 

In Fig. 2 we compare the relativistic neutrino emissivity (|66D with the non-relativistic 
energy losses, as given by formula (p^). Due to a steady decrease of the effective mass of 
the nucleon the non-relativistic formula (short-dashed curve) predicts a decreasing of the 
emissivity along with the density increase. Appearance of hyperons in the system suppresses 
the nucleon fractions and lepton abundance. Therefore at densities, where the number 
of hyperons is comparable with the number of protons, the relativistic emissivity (solid 
curve) reaches the maximum and then also has a tendency to decrease. The relativistic 
emissivity, however, is found to be substantially larger than that predicted in the non- 
relativistic approach. 

To inspect the contribution of weak magnetism and the pseudoscalar interaction we 
demonstrate two additional graphs. The long-dashed curve demonstrates the energy losses 
obtained from Eq. (|66D by formal setting $ = 0. This eliminates the pseudoscalar contri- 
bution. The dot-dashed curve is obtained by formal replacing $ = and C^j = 0, which 
eliminates both the weak magnetism and pseudoscalar contributions. A comparison of these 
curves demonstrates the weak magnetism effects. The contribution of the pseudoscalar in- 
teraction can be observed by comparing the total neutrino energy losses (solid curve) with 
the long-dashed curve, which is calculated without this contribution. We see that, due to 
weak magnetism effects, relativistic emissivities increase by approximately 40-50%, while the 
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pseudoscalar interaction only slightly suppresses the energy losses, approximately by 5%. 




Fig. 2. The left panel shows the neutrino emissivity of the direct Urea processes among 
nucleons and electrons. The right panel is the same, but for the muon reactions. The 
emissivity is shown versus the density ratio rifc/rio for the matter composition represented in 
Fig. 1. The curves begin at the threshold density. Solid curves represent the total relativistic 
emissivity, as given by Eq. (^61). The short-dashed curves are the non-relativistic emissivity 
given by Eq. (|68|) . The dot-dashed curves show the relativistic emissivity without the 
contributions of weak magnetism and pseudoscalar interaction, and the long-dashed curves 
are the emissivity without the pseudoscalar contribution. All the emissivities are given in 
units lO^^Tg erg cm~^s~^, where the temperature Tg = K. 



VII. SUMMARY AND CONCLUSION 

In the mean field approximation, we have studied the neutrino energy losses caused by 
the direct Urea processes on nucleons in the degenerate baryon matter under beta equilib- 
rium. We have shown that the direct Urea processes in a superdense matter of neutron 
star cores are kinematically allowed only due to isovector mesons, which differently interact 
with protons and neutrons. By creating the energy gap between proton and neutron spec- 
trums the isovector mesons support a time-like total momentum transfer from the nucleon, 
as required by kinematics of the reaction. In the mean field approximation, we derived 
the matrix element of the nucleon transition current, which is found to be a function of 
the space-like kinetic momentum transfer. We have calculated the neutrino energy losses 
caused by the direct Urea processes on nucleons. Our Eq. (^) for neutrino energy losses 
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exactly incorporates the effects of nucleon recoil, parity violation, weak magnetism, and 
pseudoscalar interaction. To quantify the relativistic effects we consider a self-consistent 
relativistic model, widely used in the theory of relativistic nuclear matter. The relativistic 
energy losses are up to four times larger than those given by the non-relativistic approach. 
In our analysis, we pay special attention to the effects of weak magnetism and pseudoscalar 
interaction in the neutrino energy losses. We found that, due to weak magnetism effects, 
relativistic emissivities increase by approximately 40-50%, while the pseudoscalar interac- 
tion only slightly suppresses the energy losses, approximately by 5%. The mean field Eq. 
(|5B|) may be considered as a starting point for studying of the correlation effects. 
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APPENDIX A: FIELD THEORETICAL MODELS OF NUCLEAR MATTER 

We employ a self-consistent relativistic model of nuclear matter in which baryons, B = 
n,p, S^, S+, A, interact via exchange of a, u, and p mesons [|]. In principle, the pion 
fields should be also included in the model. However, the expectation value of the pion field 
equals zero, giving no contribution to the mean fields. Therefore, only non-redundant terms 
are exhibited in the Lagrangian density 



C = J2B ^^(td^~ g^BUJ^ - ^gpBb^' ■ r) - {Mb - g^Ba) 



B 



B 

+^ (d^ad^'a - mla^) - U (a) + l{i^,d^ - rm) /, (Al) 

which includes the interaction of baryon fields B with a scalar field a, a vector field uj^ and 
an isovector field of p-meson. In the above, B are the Dirac spinor fields for baryons, 
is the isovector field of p-meson. We denote as r the isospin operator, which acts on the 
baryons of the bare mass Mb- The leptons are represented only by electrons and muons, 
/ = e~,p^, which are included in the model as noninteracting particles. The field strength 
tensors for the u and p mesons are F^^ = d^uj^, — d^u^^ and B^^, = d^hy — dyh^, respectively. 
The potential U (cr) represents the self-interactions of the scalar field and is taken to be of 
the form 

U (a) = h,M {g^^af + {g^^af . (A2) 

In what follows we consider the mean field approximation widely used in the theory of 
relativistic nuclear matter. In this approximation, the meson fields are replaced with their 
expectation values 

a^{a) = ao, u'' ^ {u^) = 0006^0, ^ {b^') = {0,0, po) 6^0- (A3) 
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In this case only the baryon fields must be quantized. This procedure yields the following 
linear Dirac equation for the baryon field 

id^Y - ^^b7°^o - ^gpBl^PoTs - (M - g^BCTo)^ B {x) = 0, (A4) 

The effective baryon mass is = Mb — QaB'^Oi and the single-particle energies and self- 
consistent potentials are given by 

Ub = 9ujB^o + gpshspo- (A5) 

Here ojq, Bq, and do are, respestively, nonzero average values of the u-, p-, and cx-meson fields, 
and gu}B, dpB, and gaB are the strong interaction couplings of the baryons to different meson 
fields, and t^B is the third component of isospin for the baryons. The number density for 
baryon species B is related to the corresponding Fermi momentum as 

- (A6) 

Thus, the average meson fields are solutions of the following equations 
^1^0 = J^dwsnB, 

B 

'nT'lpo = J29pB't3BnB, (A7) 

B 

^""^^ ^ [MbPb^Pb + - arcsm — 1 . 

For a fixed total baryon number density 

n}, = ^nB, (A8) 

B 

the additional conditions needed to obtain a solution are given by the charge neutrality 
requirement 

ne + n^ + riY,- = rip + ns+ (A9) 
and the beta equilibrium relations 

Pn = Pp + Pe, Pp = Pe, Pt,- = Pn + Pei 
PT.0 = PA= Pn, PT.+ = Pp (AlO) 

For nucleons, the strong interaction constants |^ 

^ = 3.15/^, ^ = 2.20/m, ^ = 2.19 /m, 

rria nip 

b = 0.008659, c = -0.002421, (All) 

with = 783 MeV, nip = 770 MeV, = 520 MeV, are determined by reproducing 

the nuclear matter equilibrium density no = 0.16 fni~^, and the binding energy per nucleon 
(~ 16 MeV), the symmetry energy (~ 30 — 35 MeV), the compression modulus ( 200 MeV< 
Ko < 300 MeV), and the nucleon effective mass M* = (0.6 - 0.7) x 939 MeV at no. For 
hyperons we take g^n = O.Gg^N, gpH = O.GgpN, g^H = 0.658gu,N- 
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